#model
model{
   for(i in 1:nNo){ NoImage[i] ~ dnorm(muNo, precNo)}
   for(i in 1:nYes){ Image[i] ~ dnorm(muYes, precYes)}

   sNo <- sqrt(1/precNo)
   sYes <- sqrt(1/precYes)
   diff <- muNo-muYes
   vratio <- (precNo/precYes)

   muNo~dnorm(.5, .00000000000000000001)T(0,1000000)
   muYes~dnorm(.5, .00000000000000000001)T(0,1000000)
   precNo~dgamma(.001,.001)
   precYes~dgamma(.001,.001)
}



#data
list(nNo=43, nYes=35,NoImage=c(47.2, 22.0, 20.4, 19.7, 17.4, 14.7, 13.4, 
13.0, 12.3, 12.2, 10.3, 9.7, 9.7, 9.5, 9.1, 8.9, 8.9, 8.4, 
8.1, 7.9, 7.8, 6.9, 6.3, 6.1, 5.6, 4.7, 4.7, 4.3, 4.2, 3.9, 
3.4, 3.1, 3.1, 2.7, 2.4, 2.3, 2.3, 2.1, 2.1, 2, 1.9, 1.7, 1.7
),Image=c(19.7, 16.2, 15.9, 15.4, 9.7, 8.9, 8.6, 8.6, 7.4, 
6.3, 6.1, 6, 6, 5.9, 4.9, 4.6, 3.8, 3.6, 3.5, 3.3, 3.3, 2.9, 
2.8, 2.7, 2.4, 2.3, 2, 1.8, 1.7, 1.7, 1.6, 1.4, 1.2, 1.1, 1))

#inits
list(muNo=1,muYes=1,precNo=1,precYes=1)



# Second model
#model
model {
   
   for(i in 1:nNo){ NoImage[i] ~ dgamma(shapeNo,rateNo)}
   for(i in 1:nYes){ Image[i] ~ dgamma(shapeYes,rateYes)}

   muNo <- shapeNo/rateNo
   muYes <- shapeYes/rateYes
   sdNo <- sqrt(shapeNo/rateNo)
   sdYes <- sqrt(shapeYes/rateYes)
   diff <- muYes-muNo
   vratio <- sdYes*sdYes/sdNo/sdNo

   shapeNo~dunif(0,1000)
   rateNo~dunif(0,1000)
   shapeYes~dunif(0,1000)
   rateYes~dunif(0,1000)

}

#inits
list(shapeNo=1,rateNo=1,shapeYes=1,rateYes=1)